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The recent WMAP data have confirmed that exotic dark matter together with the vacuum energy 
(cosmological constant) dominate in the flat Universe. Supersymmetry provides a natural dark 
matter candidate, the lightest supersymmetric particle (LSP). Thus the direct dark matter detection 
is central to particle physics and cosmology. Most of the research on this issue has hitherto focused 
on the detection of the recoiling nucleus. In this paper we study transitions to the excited states, 
' focusing on the first excited state at 50 keV of Iodine A=127. We find that the transition rate to 

, this excited state is ^ 10 percent of the transition to the ground state. So, in principle, the extra 

f^**) . signature of the gamma ray following its de-excitation can be exploited experimentally. 

<N , 

^ PACS numbers: 95.35.+d, 12.60.Jv 11.30Pb 21.60-n 21.60 Cs 21.60 Ev 

O ' 

O 

oo 

The combined MAXIMA-1 |l|, BOOMERANG |2j, DASI [3|, COBE/DMR Cosmic Microwave Background (CMB) 
observations 0, the recent WMAP data [|J and SDSS imply that the Universe is flat |j| and and that most of the 
^ . matter in the Universe is dark, i.e. exotic. 

in 

CO 

for baryonic matter , cold dark matter and dark energy respectively. An analysis of a combination of SDSS and 
! WMAP data yields n m « 0.30 ± 0.04(1<t). Crudely speaking 

o ■ 

^1 6 « 0.05, Qcdm « 0.30, Oa « 0-65 

^ Since the non exotic component cannot exceed 40% of the CDM ||, there is room for the exotic WIMP's (Weakly 

Q-i, Interacting Massive Particles). In fact the DAMA experiment has claimed the observation of one signal in direct 
detection of a WIMP, which with better statistics has subsequently been interpreted as a modulation signal [ic| . 
although these data are not consistent with other recent experiments, see e.g. EDELWEISS |E3 and CDMS 0. 

Supersymmetry naturally provides candidates for the dark matter constituents p^ | .[l ^| -|l7 j . In the most favored 
scenario of supersymmetry the LSP can be simply described as a Majorana fermion, a linear combination of the 
$_i ' neutral components of the gauginos and higgsinos 0,0]-[2l]. In most calculations the neutralino is assumed to be 
primarily a gaugino, usually a bino. Models which predict a substantial fraction of higgsino lead to a relatively large 
spin induced cross section due to the Z-exchange. Such models tend to violate the LSP relic abundance constraint 
and are not favored. Some claims have recently been made, however, to the effect that the WMAP relic abundance 
constraint can be satisfied in the hyperbolic branch of the allowed SUSY parameter space, even though the neutralino 
is then primarily a higgsino j^- We will not elaborate further on this point, but we will take the optimistic view that 
the detection rates due to the spin may be large enough to be exploited by the experiments, see, e.g., |24|. 
|25| . Such a view is further encouraged by the fact that, unlike the scalar interaction, the axial current allows one 
to populate excited nuclear states, provided that their energies are sufficiently low so that they are accessible by the 
low energy LSP, a prospect proposed long time ago by Ejiri and collaborators j2||. As a matter of fact the average 
kinetic energy of the LSP is: 



n b = 0.044 ± 0.04, fi TO = 0.27 ± 0.04, n A = 0.69 ± 0.08 



So for sufficiently heavy LSP the average energy may exceed the excitation energy, e.g. of about 50 keV for the 
excited state of 127 1. In other words one can explore the high velocity window, up to the escape velocity of v esc « 
570 km/ s. From a Nuclear Physics point of view this transition is not expected to be suppressed, since it is of the 
type (5/2)+ (7/2)+, i.e. Gamow-Teller like. 
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THE LSP-NUCLEUS CROSS SECTION 

The LSP-nucleus differential cross section with respect to the energy transfer Q for a given LSP velocity v due to 
the spin can be cast in the form 



du 



da(u,v) = nf ^ A0 Y> spin Fii(u) (2) 



where F\\ is the spin response of the isovector channel |27j. We have used a dimensionless variable u, proportional to 
Q, which has been found convenient for handling the nuclear form factor pjll ] F(u), namely: 

u = 7T > = . 1 , 2 « 40 x A" 4 / 3 MeV. (3) 
Qo Am N b z 

fi r is the reduced LSP-nucleus mass and b is the (harmonic oscillator) nuclear size parameter. 
The quantity £ sp m is given by: 

^spin ( 7 \") ^-n ^0 CspirojCspm To ^(^0 (4) 

Mr(p) P,X 3(1 + 4)2 

J A 

firiv) ~ mp is the LSP-nucleon reduced mass and cr^™ is the proton cross-section associated with the spin and S(u) 
|27j is given by: 

= [(^o)) 2 ^ + 2 ^no(o)nx(o)^ + fMo)) 2 ] (5) 

t\x \u) J A txi{U) 



The overall normalization of the axial current components is not important in the present discussion. We mention, 
however, that they have been normalized so that tr^" = 3(/^ + F\) 2 a , c* p *o = 3(/^ — F\) 2 ao, cr = ^(G^ m p ) 2 , 
for the proton and neutron spin cross sections respectively. 

Notice that S(u), being dependent on ratios, is expected to be less dependent on the energy transfer and the 
scale of supersymmetry. As a matter of fact it has been found [13 that the spin response functions Fij , the " spin 
form factors" normalized to unity at u = 0, have similar dependence on the energy transfer. In other words S(u) is 
essentially independent of u. 

A number of nuclear spin matrix elements for a variety of targets have become available |27l | -|36j | . Some static spin 
matrix elements [2]), [23, for some nuclei of interest are given in table [I] The spin matrix elements appearing in 
the table are defined as: 

fio = S p + £„ ,Qi = £ p -£ n , E fc = —====< J/\\2 Sk\\Ji >, k=p,n, (6) 

V 2Ji + 1 

where by double bar we indicate the reduced matrix element as defined in standard textbook as e.g the one by 
Edmonds. The < Sk > , fc = p,n, are defined in the literature, e.g. Ressel and Dean |2^|. 

In order to proceed, in addition to the static spin matrix elements and the spin response functions Fqq, Fqi, Fh 
("spin form factors"), normalized to unity at zero momentum transfer, one must know the isoscalar and the isovector 
axial currents. In the allowed supersymmetric parameter space one obtains the isoscalar and the isovector axial 
current components at the quark level. Going from the quark to the nucleon level is straightforward for the isovector 
current, but it not trivial in the case of the isoscalar current, since the naive quark model fails badly. Most of the 
proton spin is not due to the quark spins (proton spin crisis-EMC effect). Thus one finds: 

f%(Q)->f% = 9°Afl(Q) >f\(Q)-> & = >9°a « 0.1 , 3^ = 1.23 (7) 

It thus appears that, barring very special circumstances, whereby the isoscalar contribution at the quark level, /5(<z), 
is much larger than the isovector, f\{q), the isoscalar contribution can be neglected. 

THE STRUCTURE OF 127 1 



This nucleus has a ground state 5/2 + and a first excited state a 7/2 + at 57.6keV. As it has already been mentioned 
it is a popular target for dark matter detection. As a result the structure of its ground state has been studied 
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theoretically by a lot of groups. Among them we mention again the work of Ressel and Dean 29], the work of Engel, 
Pittel and Vogel [35]], Engel and Vogel (3^], Iachello, Krauss and Maiano |3^|, Nikolaev and Klapdor-Kleingrothaus 
[3^ and more recently by Suhonen and collaborators (3(|. In all these calculations it appears that the spin matrix 
element is dominated by its proton component, which in our notation implies that the isoscalar and the isovector 
components are the same. In these calculations there appears to be a spread in the spin matrix elements ranging from 
0.07 up to 0.354, in the notation of Ressel and Dean |29j. This, of course, implies discrepancies of about a factor of 
25 in the event rates. 

In the present work we are primarily interested in the spin matrix element connecting the ground state to the first 
excited state. As we have mentioned, however, it is advantageous to compute the branching ratio. In addition to 
factoring out most of the uncertainties connected with the SUSY parameters and the structure of the nucleon, we 
expect the ratio of the two spin matrix elements to be more reliable than their absolute values. Thus this ratio is 
going to be the most important factor in designing the experiments to detect the rate to the excited state. 

In order to have for Z = 53 single particle orbits suitable for the above states, close to the Fermi surface [22|, namely 
[402,5/2] and [404,7/2] proton orbits, one should consider oblate deformation j^. The value of the deformation, 
however, is not very well defined. To this end one should explore the spectrum of this system and combine information 
obtained from it with data on the neighbor ing, systems Te and Xe. On notices that on the above band heads A J = 2 
rotation aligned coupling bands are built [39j. From the spectra one may also conclude that the 7/2 + data most 
likely correspond to a deformed band. The 5/2 data correspond to a highly perturbed band or to a spherical vibrator. 
Other people have considered triaxial shapes, see ref. 9 of Shroy et al 39]. So we will take it as an educated guess 
that a description in terms of a rotor plus a particle with relatively large oblate deformation provides a reasonable 
description of this nucleus. 

Given the above framework one may eventually have to do a self-consistent calculation, which will yield both the 
single particle states and the equilibrium deformation. At this exploratory stage, however, we were merely interested 
in testing the rather crude collective model contentions discussed above. To that effect we found it appropriate to 
take as simple, and easy to read, as possible the single particle substrate of our description. Thus to estimate the spin 
matrix elements, in the present simple spectroscopic study, we have proceeded in the following oversimplified fashion: 

• Make a big leap in history to Nilsson's thesis [s^R ] . 

• Compute the spin proton matrix elements involving the above bare single particle states for three choices of 
oblate deformations. 

• include the effect of the Coriolis coupling for the 7/2 + . 
The relevant deformation parameter 77 is defined |37l ] as 



where 8 is the usual Nilsson deformation parameter and x ~ 0.05. The values of rj considered are —6, —4, —2. To first 
order the corresponding 8 values are: —0.3, —0.2, —0.1. The corresponding states in the standard notation N£AT, are 
given by: 




(8) 



X 



|5/2> 



a|442+ > +6|422+ > +c|443- > 



(9) 



V»+6 2 + c 2 



|7/2> 



d|443+ > +e|444- > 



(10) 



Vd 2 + e 2 



The numerical values for these coefficients in the above order of the 77' s are: 



rj = -6 -> (a, b, c, d, e) = (1.000, 0.552, -0.809, 1.000, -1.192) 



77 = -4 -> (a, b, c, d, e) = (1.000, 0.558, -1.081, 1.000, -1.662) 



77 = -2 -> (a, b, c, d, e) = (1.000, 0.517, -1.451, 1.000, -2.219) 
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Thus in the standard angular momentum notation \J,M > one obtains: 

< 5/2, 5/2|s z |5/2, 5/2 >= 0.166, 0.033, -0.124 (11) 

< 7/2, 7/2\s z \7/2, 7/2 >= -0.087, -0.234, -0.331 (12) 

< 5/2, 5/2|s+|(7/2, 5/2) >= 0.171, -0.120, -0.073 (13) 

The array of numbers in the above three equations corresponds to r\ = —6, —4, —2 respectively. 

One may also have to consider the effect of the Coriolis coupling, see e.g. |4jj, between the I — 7/2 and the 
K = 5/2, 7/2 bare single particle states with the Hamiltonian matrix as follows: 

Ma = E\ are + — [1(1 + 1) - 2K 2 ] (14) 



for the diagonal term. For the off-diagonal term one has: 

V = 



h 2 r/(/ + l) -K < (K < + 1) 



J 



1/2 

(X>|j+|X<) (15) 



where j + is given in units of h . One, of course, has to adopt a reasonable value for the moment of inertia J. For 
such a choice and rj — —6 one finds: 

Mn = 387, Af 22 = 145, V = 207 (16) 

(in KeV). The resulting components of the eigenfunctions are: 

0.50,0.87 and -0.87,0.50 for the 5/2 and 7/2 respectively. 
With the above wave functions we find the ground state spin matrix elements: 

< S p >= 0.166 , < S„ >= 

, which are in good agreement with Iachello et al [3^ and Nikolaev and Klapdor-Kleingrothaus [34| , but substantially 
smaller than Ressel and Dean |2^|. In our notation this leads to: 

fll = n\ = fio^i = 0.164 (17) 

This value is also smaller than the more recent result 36] . For the transition to the excited state we have obtained: 

fl 2 = fl 2 = fi f7i = 0.312 (18) 

Thus the spin matrix clement, left to itself, seems to yield an enhancement of about a factor of two for the transition 
to the excited state. In our simple model this result is independent of the ratio of the isoscalar to isovector amplitude. 

EXPRESSIONS FOR THE RATES 

In this section we will briefly outline the expressions yielding the event rates, both modulated and unmodulated, 
(for more details see our earlier work j2j| and references therein). The differential non directional rate can be written 
as 

dR= P(0) A m da(u,v)\v\ (19) 
m x ArriN 

where da(u, v) was given above, p(0) is the LSP density in our vicinity, m x the LSP mass, m is the mass of the target 
and A is the nuclear mass number. 

For a given velocity distribution f(t/), with respect to the center of the galaxy, one can find the velocity distribution 
in the lab frame f(v, VE) by writing 

v'=v+ V E , V E =V + Vl 
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where Vo = 229 km/ s is the sun's velocity (around the center of the galaxy) and V\ the Earth's velocity (around 
the sun) . The velocity of the earth is given by 

Ve = vqz + i>i( sina x — cosa cosj y + cosa sirvy z ) (20) 

In the above formula z is in the direction of the sun's motion, x is in the radial direction out of the galaxy, y is 
perpendicular in the plane of the galaxy (jj — z x x) and 7 ss 7r/6 is the inclination of the axis of the ecliptic with 
respect to the plane of the galaxy, a is the phase of the Earth in its motion around the sun (a = around June 2nd) . 
The time dependence of the event rate, modulation effect, will be present for transitions to the excited states as well. 
Folding with the velocity distribution we get: 

dRundir \ _ /dR\ _ p(0) m /r^r.g. 
du I \ du I m x ArriN du 



where 



<£-/■ 



rLffrvJ^^fv (22) 
In our calculation we have used the standard M-B distribution: 

/<»') = j^m^-^f) m 

In what follows we will express the LSP velocity in units of vq introducing the dimensionless quantity y — v/vq. 
Incorporating the relevant kinematics and integrating the above expression, Eq. 1221) . from u m i n to u max we obtain 
the total rates as follows: 

R °- = £."^ du (24) 

with u min = Q m i n /Qo, Q m in imposed by the detector energy cutoff and u max = (y esc /a) 2 , with a = [^ r frV2]~\ is 
imposed by the escape velocity (y esc — 2.84). 

Rexc = I 11 ™* (^(1 " (25) 

Ju^ a \ du u z I 

where u exc — / LrE % l and E x is the excitation anergy of the final nucleus. The upper limit is now given by u max — 
{y/af - (E x /Q ) 

BRANCHING RATIOS 

In the present calculation, to spare the reader with details about the SUSY allowed parameter space, we are not 
going to discuss the absolute value of the event rates. We will, instead, estimate the ratio of the rate to the excited 
state divided by that to the ground state (branching ratio of the rates) as a function of the LSP mass. This is of 
interest to the experimentalists. It can be cast in the form: 

J3RR ^e%c(fy ^ exci.'U'exci Uumax)^- ~t~ h eX c{v>exci'U i max) COS Ql\ (26) 

S gs (0) ^ g S {u mm )[l + h(u min ) cosa] 

In the above expression h and h exc are the modulation amplitudes for the ground state and excited state respectively. 
They correspond to the ratio of the amplitude associated with the motion of the Earth divided by the one which is 
independent of the motion of the Earth. The parameter a is the phase of the Earth in its motion around the Sun 
(a = on June 3nd) . The modulation amplitudes can only be obtained numerically and are not going to be discussed 
in detail here (for more information see our earlier work |25J). 

S gs (0) and S exc (0) are the static spin matrix elements evaluated via Eq. [S] This ratio of the static spin matrix 
elements is essentially independent of supersymmetry, if the isoscalar contribution is neglected due to the EMC effect. 
The functions \? are given as follows : 

* gs (u min )= ^g<(«)[^-4. e )] du (27) 
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FIG. 1: The ratio of the rate to the excited state divided by that of the ground state as a function of the LSP mass (in GeV) 
for 127 1. We assumed that the static spin matrix element of the transition from the ground to the excited state is a factor of 
1.9 larger than that involving the ground state, but the spin response functions Fn(u) are the same. On the left we show the 
results for Qmin = and on the right for Qmin = 10 KeV. 



$e IC (« ac , u max ) = [ ?pM F g c (u)(l - ^f) ty(aVS(l + W«)) - lKife«0] du (28) 

The functions ip were evaluated numerically, with the correct kinematics imposed on the neutralino velocity in the 
laboratory, which is obtained by translating the simple relation v < v esc , which holds with respect to the galactic 
center. To a very good approximation, however, they can be given by: 

iP(z) = ^[Erf(z+l)-Erf(z-l)] (29) 

and with Erf(x) the well known error function: 

Erf{x) = —= / e~* 2 dt. 
V 71 " Jo 

As we have mentioned the dependence of S(u) on the energy transfer is extremely mild and it can be neglected. 



CONCLUSIONS 

In the present paper we have calculated the branching ratio BBR for the neutralino - 127 I scattering event rate to 
the first excited state. We found that the ratio of the static spin matrix element involving the excited state is 1.9 times 
larger than that of the ground state. Assuming that the spin response functions Fu(u) are identical, we obtained 
the ratio BRR (branching ratio of the rate to the excited divide by that to the ground state), which is exhibited in 
Fig. ^ In computing these branching ratios we did not include the quenching factors affecting the nuclear recoils, 
i.e. the detected event rates for transitions to the ground state. If such factors are included, the branching ratios will 
increase. 

The difference in the branching ratios of Fig. ^ can be understood as follows: Since the transitions the excited 
state will be detected not by nuclear recoil, but via the EM decay of the final state, they do not suffer from the 
threshold energy cut off. This cut off suppresses only the transitions to the ground state, which are detected through 
nuclear recoils. Anyway the branching ratio is affected. Thus this ratio is very small in the, at present, unattainable 
experimentally case of Q m in — 0. As we have mentioned, however, even in this case, the branching ratio will increase, 
if the quenching factor is taken into account. The branching ratio is small even in the case of Q m in = 10 keV for LSP 
masses less than 50 GeV. This is not unexpected, since such a light neutralino cannot provide the energy needed to 
excite the nucleus. For Q m in — 10 this branching ratio attains a maximum value, around 10%, for LSP masses around 
100 GeV and gets values around 6% for larger masses. For larger Q m in the branching ratio will further increase. 

The smallness of the branching ratio, < 10%, should not discourage the experiments, since the transition to the 
excited state offers experimental advantages. This is especially true, if the traditional experiments cannot reach 
Qmin d 10 keV. 

Finally from Figs [5] and [3] we notice that, especially for m x > 100 GeV, the modulation amplitude for transitions 
to the excited state is much bigger compared to that of the elastic scattering. We hope that this additional signal 
can perhaps be exploited by the experimentalists to discriminate against background 7 rays, in a fashion analogous 
to the ongoing measurements of nuclear recoils. 
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FIG. 2: The the same as in Fig. Qfor the modulation amplitude h for the transition to the ground state. For the definitions 
see text. 




FIG. 3: The the same as in Fig. |2]for the modulation amplitude h exc for the transition to the excited state. 
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